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We report on the successful demonstration of selective acceleration of deuterium ions by target-normal sheath 
acceleration (TNSA) with a high-energy petawatt laser. TNSA typically produces a multi-species ion beam 
that originates from the intrinsic hydrocarbon and water vapor contaminants on the target surface. Using 
the method first developed by Morrison, et al.,^ an ion beam with >99% deuterium ions and peak energy 
14 MeV/nucleon is produced with a 200 J, 700fs, > lO^^W/cm^ laser pulse by cryogenically freezing heavy 
water (D 2 O) vapor onto the rear surface of the target prior to the shot. Within the range of our detectors 
(0-8.5°), we find laser-to-deuterium-ion energy conversion efficiency of 4.3% above 0.7 MeV/nucleon while a 
conservative estimate of the total beam gives a conversion efficiency of 9.4%. 


I. INTRODUCTION 

The acceleration of ions with high-power lasers has 
drawn a great deal of attention over the last 15 years. 
This interest has been driven by a wide range of promis¬ 
ing applications^ coupled with the constantly improving 
capabilities of laser facilities. Bunches of MeV ions can 
be used for radiobiological studies possibly relevant for 
cancer therapy,^’^ creation of warm dense matter with 
isochoric heating,^ proton driven fast ignition,^ and neu¬ 
tron production. The best known and most investi¬ 
gated mechanism for generating fast ion beams with a 
high-intensity laser is target-normal sheath acceleration 
(TNSA),^’^ which produces protons with lO’s MeV ener¬ 
gies. The ions are accelerated by the sheath field that 
is formed on the rear surface of the target by laser¬ 
generated hot electrons.^ TNSA ion beams typically 
have a broad distribution in energy and emission an¬ 
gle with laser-to-ion conversion efficiency up to several 
percent. 

Of particular recent interest is the generation of neu¬ 
tron beams from laser-irradiated targets; one mechanism 
for laser-based neutron production is the pitcher-catcher 
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scheme.Here, the laser produces an ion beam (typically 
protons or deuterium ions) from a primary “pitcher” tar¬ 
get. The ions are collided into a secondary “catcher” 
target that is composed of a fusable material. There are 
several reactions which can produce neutrons in this con¬ 
figuration including p(Li,n), d(d,n), d(t,n), d(Li,n), and 
d(Be,n). The deuterium cross sections are all larger than 
p(Li,n) but, without special care, TNSA predominantly 
accelerates protons; until recently, progress in develop¬ 
ing a deuterium source that can utilize the larger cross 
sections has been limited. 

There are several emerging, alternative mechanisms 
for the production of intense particle beams from laser- 
plasma interactions: e.g., radiation pressure acceleration 
(RPA)^^’^^ and breakout afterburner (BOA).^^ However, 
these mechanisms have strict laser and target require¬ 
ments and are still under active development. Currently, 
the simplest configuration to produce a deuterium beam 
for neutron generation is TNSA from a deuterium-rich 
target (deuterated plastic, for example). However, the 
ion beams produced by this configuration are ubiqui¬ 
tously dominated by C+, 0+, and apparently, 

controlling the accelerated ion species takes more than 
target material selection. The reason for this is well 
known: sub-micron layers of hydrocarbon and water- 
vapor contaminants cover the targets. The difficulty in 
satisfactorily employing the pitcher-catcher technique for 



2 


neutron production has been discussed by Willingale, et 
al.^^ 

There have been several attempts at reducing the yield 
of contaminant ions in TNSA. Heating a 50 /im thick A1 
target to 600 K has been shown to reduce the peak en¬ 
ergy and yield of protons by about an order of magnitude 
while enhancing species coated onto the substrate. 
Ablation with a secondary laser^^ has also been shown to 
reduce the proton signal by decreasing the rear surface 
field. Unfortunately, none of these contaminant reduc¬ 
tion schemes have produced the desired quality of deu¬ 
terium ion sources. 

Here we report the results of an experiment that ex¬ 
tends the method developed by Morrison, et al.,^ to a 
new regime of laser energy and intensity. This approach 
freezes a /im’s thick layer of heavy water over the ubiqui¬ 
tous proton-rich contaminants. Morrison demonstrated 
that this method produces an ion beam with ^ 99% deu¬ 
terium ions while maintaining typical laser-to-ion conver¬ 
sion efficiency. The ice layer addition can be done quickly 
and reproducibly, and is synchronized with the laser to 
prevent regrowth of contaminants. We have produced an 
ion beam that has a similar deuterium-to-proton ratio to 
the heavy water (>0.99) and high laser-to-deuterium-ion 
conversion efficiency. 


II. EXPERIMENTAL SETUP 

This experiment was performed using the petawatt 
arm of the Vulcan laser at the Rutherford Appleton Lab¬ 
oratory (RAL). As configured for this experiment, the 
laser delivered 200 J on target in a 700 fs pulse of 1053 
nm light and is focused to above 10^^ W/cm^ with an //3 
off-axis parabola to a ^6 fim full-width at half-maximum 
spot. The laser is normally incident off of a plasma mir¬ 
ror onto the front surface of a 10 jam thick An foil with 
a 3 fim (±1 /im) thick layer of frozen heavy water on the 
rear surface. 


A. Ice Formation 

Figure 1 shows the in-chamber portion of the setup 
used in this experiment. The ice layer is formed by cryo- 
genically cooling the An target and releasing a small puff 
of heavy water vapor just before the laser shot. The va¬ 
por source is a trapped volume of heavy water that had a 
measured ambient vapor pressure of 11 Torr. The outlet 
nozzle was positioned ^6 cm from the target at a ^50° 
angle with respect to the target surface. The vapor is let 
into the chamber just before the shot with a solenoidal 
valve that is synchronized with the the laser system. 
For the data presented, we used the maximum possible 
growth rate of ^ l/rm/s, found empirically to produce 
the best deuterium beams. ^ 3/im was the minimum 
achieveable thickness due to technical reasons related to 
triggering. Thinner ice layers, which may produce even 


higher TNSA laser-to-ion conversion efficiency,^^’^^ could 
be readily produced with improved timing. 


I = D 2 O ice {3+/- liim thick) 

I = Cold Au target (lO^m thick) 
= Cold target post 
^ = Plastic insulating coupler 
I = Target wheel 

Not drawn to scale 


Figure 1. (color online). The target configuration for freez¬ 
ing heavy water on a 10 /xm thick Au foil target. The laser 
pulse reflects off of the plasma mirror and is normally inci¬ 
dent on the target (shown in yellow). The heavy water vapor 
(blue) comes in through a standard 6 mm outer diameter tube 
(black) whose outlet is ~6 cm from the Au target. The target 
is mounted with an A1 post that is thermally isolated from 
the vacuum chamber and target stage (green) by an insulat¬ 
ing coupler (black). The target is cooled to below -100 C by a 
liquid nitrogen reservoir (purple) which is flexibly connected 
(also purple) to the target post. 

The heavy ice thickness was characterized with dy¬ 
namic thin-film interference reflectometry. The puff 
duration-to-ice thickness calibration was done pre-shot 
- a conflict of diagnostics within the target chamber pre¬ 
vented in situ characterization. A visible cw laser is re¬ 
flected off of the surface and onto a charge-coupled de¬ 
vice. As the ice grows, constructive or destructive in¬ 
terference peaks occur when the optical path is equal 
to alternating half-integer multiples of the laser wave¬ 
length. Counting the number of these peaks that occur 
during the growth gives the total thickness. Critically, 
the targets are coated with a surfactant which enables 
optically smooth ice formation. The surfactant is applied 
outside of the target chamber before the shot and opti¬ 
cally smooth ice can be grown even after several hours 
in vacuum. Thorough pre-shot testing indicated that the 
ice layering was repeatable for the same conditions of 
heavy water vapor pressure in our apparatus and growth 
duration. 

B. Ion Spectra Characterization 

The ion spectra are recorded with 4 Thomson parabola 
spectrometers (TPS) at different angles (-6°, 0°, 3.5°, 
and 8.5°); there was an additional TPS at 30° that 
recorded no ion signal. Each TPS has a 100 /xm radius 
pinhole which subtend between 1.8-2.1x10“^ sr. The 
ions are detected with BAS-TR^^ image plates (IP). A 
previously ubiquitous problem in this type of experiment 
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a) TPS layout 4X TPS 



b) Sample IP filtering 



Increasing energy 


Figure 2. (color online). The detector arrangement is shown 
in (a); there was an additional TPS at 30° that recorded 
no ion signal. The 100 fim. radius pinholes for the TPS are 
all located between 1.23-1.31 m from the target and subtend 
between 1.8-2.1x10“^ sr. An example of the differential IP 
filtering is shown in (b). The filters are chosen such that 
heavier ions like and are stopped before reaching the 
IP while light ions (H^ and D^) are tramsitted with reduced 
kinetic energy. The entire IP is covered in 6 /xm thick A1 so 
the region labeled 12 /xm A1 is hltered with a total of 18 /xm 
Al, etc. 


is segregating the proton and deuterium ion signals from 
the higher mass contaminants. There are two causes for 
this. First, the TPS disperses ions by their charge-to- 
mass ratio and so, for example, the deuterium ion (D^+), 
and etc. ion signals are overlapped. Sec¬ 

ond, the highest energies of neighboring charge-to-mass 
species can be overlapped due to insufficient dispersion 
and non-zero pinhole diameter. Here, the overlapping 
of heavier ions with the deuterium and proton signals is 
prevented by using differential filters for the IPs. 

The layout and an example of the differential filtering 
is shown in Figure 2. The details of the filtering are dis¬ 
cussed thoroughly by Alejo, et al.^^ The thicknesses and 
materials of the IP liters are chosen to stop the higher- 
mass contaminant ions while allowing the deuterium ions 
and protons through to the IP. The entire IP is covered 
with a 6 jam thick Al foil so the 12 jam Al filter region is 
actually filtered by 18 jam thick Al and so on. 

The energy dispersion onto the image plates is cal¬ 
culated numerically using a previously benchmarked 
method^^ which accounts for fringe fields that are nor¬ 
mally not taken into account for this type of diagnostic. 
The TPS consists of two regions of separate and parallel 
static magnetic and electric fields which are achieved by 
a yoked permanent magnet pair (peak 1.02 T) and paral¬ 
lel conducting plates with a large potential difference (16 


kV, 15mm separation). We use RADIA,^^ a three dimen¬ 
sional magnetostatic solver, to characterize the magnetic 
fields - basic features are double checked with a Teslame- 
ter. Laplace’s equation is solved using Matlab’s built in 
partial differential equation solver which gives the elec¬ 
tric field in two dimensions. Test particles are numer¬ 
ically propagated through this system and are spatially 
mapped onto the IP giving the energy dispersion for each 
TPS. The absolute response and decay properties of the 
IP are described in Alejo, et al.^^. 

In some cases, the deuterium signal was so strong that 
the scanner recorded regions of saturation. The PSL val¬ 
ues for these pixels are calculated by rescanning the IP 
to generate a calibration function. The saturated values 
are extrapolated using this function; an example of this 
is in Figure 3, which shows the correction for the scan 
shown below in Figure 4. 



Figure 3. (color online). Sample correction curve for a satu¬ 
rated scan: the unsaturated pixels are shown as blue dots and 
the saturated pixels are shown as red diamonds. The PSL val¬ 
ues for unsaturated pixels are used to construct a quadratic 
fit which is extrapolated to determine the actual signal. The 
correlation coefficient (R^) for the ht is 0.98. 


III. RESULTS AND DISCUSSION 

Raw PSL counts (logio scale) from a sample shot (0° 
spectrometer) is shown in figure 4. The bright parabolic 
line is the deuterium signal (labeled D+); below it a very 
faint proton line (labeled H+) is visible. The filters (Cu, 
Fe, Al 70 /xm, etc.) used in each region are labeled at the 
top. Between each filter there is a gap of approximately 
1mm; in the Fe/Al 70/xm and Al 70/xm/Al 25 jam gaps 
there is a faint 0+ track visible. The three other TPS 
recorded quantitatively similar results; strong deuterium 
ion signal and nearly extinguished contaminant signal. 

The absolute energy spectra from all 4 detectors are 
shown in Figure 5 (same shot as Figure 4). The maxi¬ 
mum energy deuterium ion, recorded by the 8.5° TPS, 
was 14 MeV/nucleon. The deuterium ion-to-proton ratio 
for the spectrometer in Figure 4 is better than 0.99, con¬ 
sistent with the purity of the source heavy water. Across 
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Raw PSL 0 Thomson (log^^ scale) 



Figure 4. (color online). Raw PSL counts from a sample shot. The bright yellow quadratic track is the deuterium ion signal 
(labeled D+). The extremely faint proton signal (labeled H+) is below the deuterium track. The regions corresponding to the 
differential Liters are labeled (Cu, Fe, etc.) at the top. There are gaps between each Liter - faint tracks are visible in the 
Fe/Al 70/im and the A1 TO/xm/Al 25/im gaps. 


5 shots (20 recorded spectra), a ratio of ^0.9 or better 
is typical and the minimum observed ratio is ^0.7. The 
differential Liters blind the detectors to heavier ions but 
one can assume that these heavier ions are also mini¬ 
mized. This assumption is based on the relatively low 
signal observed in the Liter gaps and the minimal proton 
signal. Protons should be preferentially accelerated over 
other species because they are the lightest ions of all can¬ 
didates; their absence is strong evidence for minimization 
of other contaminant species. 


Deuterium ion beam energy spectrum 



Figure 5. (color online). Energy spectra recorded by the 
TPS positioned at 8.5° (blue), 3.5° (yellow), 0° (red), and 
-6° (green). The horizontal axis is given in MeV/nucelon. 

Finally, we estimate the deuterium ion beam properties 
above 0.7 MeV/nucleon, the minimum energy incident on 
the detector. Inside our detector limit, the conversion ef- 
Lciency is 4.3%, found by assuming azimuthal symmetry 
and making a linear Lt to the TPS over the polar angle. 


However, this is likely a signiLcant underestimate due to 
the angular cut-off. The full-beam conversion efficiency 
can be estimated by considering the typical divergence 
of ions reported in literature for similar conditions. For 
instance, Maksimchuk et ah, using a similar heavy ice 
target conffguration with a much smaller laser (6 J, 15 
TW), observed a Gaussian angular distribution with a 
10° Using these results, we Lnd a conversion 

efficiency of 9.4%. On the other hand, Nurnberg et al.^^ 
observed a proton beam that was much broader than 
Maksimchuk. These results used the same plasma mirror 
conffguration at Vulcan PW as ours suggesting that 9.4% 
efficiency could be an underestimate. We are planning to 
characterize the full energy-resolved angular distribution 
as well as repeatability in a future experiment. 


IV. CONCLUSIONS 

We have demonstrated the ability to produce a nearly 
pure deuterium ion beam from cryogenic An targets 
coated by a layer of heavy ice with a high-energy petawatt 
laser. Inside our detector limit (0 - 8.5°) we observe 
an ion beam with >0.99 deuterium-to-proton yield ratio, 
high peak energy (14 MeV/nucleon), and high conversion 
efficiency (4.3%); a conservative estimate for the total 
conversion efficiency is 9.4%. Further investigations will 
be done to fully characterize the deuterium ion beam in 
the future. 
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